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ABSTRACT 

We investigate the evolution of molecular abundance in quiescent protoplanetary 
disks which are presumed to be around weak-line T Tauri stars. In the region of surface 
density less than 10 2 g cm -2 (distance from the star ^ 10 AU in the minimum-mass 
solar nebula), cosmic rays are barely attenuated even in the midplane of the disk and 
produce chemically active ions such as He" 1 " and Hjj". Through reactions with these ions 
CO and N2 are finally transformed into CO2, NH3, and HCN. In the region where the 
temperature is low enough for these products to freeze onto grains, considerable amount 
of carbon and nitrogen is locked up in the ice mantle and is depleted from the gas phase 
in a time scale ^ 3 x 10 6 yr. Oxidized (C 0 2 ) ice and reduced (NH 3 and hydrocarbon) 
ice naturally coexist in this part of the disk. The molecular abundance both in the gas 
phase and in ice mantle varies significantly with the distance from the central star. 


Subject headings: circumstcllar matter — ISM: molecules — stars: formation — stars: 
pre-main sequence — solar system: formation 
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1. INTRODUCTION 

Radio, infrared, and optical observations have re¬ 
cently revealed the existence of circumstellar disks 
around young stellar objects (Sargent & Beckwith 
1994; O’dell & Wen Zheng 1994). Planet formation 
must be going on at least in some of them. Theoreti¬ 
cal study on the evolution of molecular abundance in 
such “ protoplanetary ” disks is important from var¬ 
ious points of view. First, it will directly show what 
material the bodies in the planetary systems are made 
from. Secondly, molecular abundance can be a use¬ 
ful probe in investigating the formation processes of 
planetary systems. For example, we may be able to 
reduce the processes, places, and epochs of the forma¬ 
tion of primitive bodies, such as comets, by comparing 
their molecular composition with the theoretical re¬ 
sults on the distribution and evolution of molecular 
abundance in the disk. Thirdly, study of the structure 
and evolution of the disk by observations of molecu¬ 
lar lines should ask for a help of theoretical study on 
molecular evolution. For example, the amount of the 
gaseous component in the disk is estimated by obser¬ 
vations of molecules other than the main component 
H 2 because H 2 can hardly emit photons except at very 
inner region of the disk. Knowledge on the abundance 
of the molecules relative to hydrogen is indispensable 
for such study. 

Some theoretical works have been done on the 
molecular abundance in protoplanetary disks. Prinn 
and his colleagues (Prinn 1993, and references therein) 
investigated the molecular evolution in fully turbu¬ 
lent accretion disks and proposed the kinetic inhibi¬ 
tion (KI) model. Their basic idea is that the molec¬ 
ular abundance in the matter flowing outwards is “ 
quenched ” when its temperature has decreased to a 
value below which the time scale for chemical reaction 
is larger than the dynamical time scale. Adopting the 
lifetime of the disk ~ 10 13 s, they concluded that the 
molecular abundance was quenched at 850 — 1500K, 
and that CO and N 2 were the dominant component 
in the solar nebula. They considered neutral-neutral 
reactions which can be efficient in the region of T 
several xlO 2 K. 

In this letter we investigate the evolution of molec¬ 
ular abundance in a relatively quiescent disk in the 
post-accretion phase, the so-called weak-line T Tauri 
(WTTS) phase and thereafter, in which transport and 
mixing of matter is not so efficient as in the earlier 
active phases mentioned above. We consider the disk 


regions where the column density is lower than the 
attenuation length of cosmic rays, 96 g cm -2 (Ume- 
bayashi & Nakano 1981). The minimum-mass solar 
nebula (Hayashi 1981) satisfies this condition in the 
region more than about 7AU away from the central 
star with the temperature T <; 10 2 K. We take into 
account the ionization by cosmic rays and the subse¬ 
quent ion-molecule reactions, which were not consid¬ 
ered in the previous works. 

2. DISK MODEL AND REACTION NET¬ 
WORK 

As the model of the protoplanetary disk, we adopt 
the so-called minimum-mass solar nebula (Hayashi 
1981), which has the the distribution of the surface 
density T,(R) = 54 (i?/10AU)^ 3 ^ 2 g cm -2 and tem¬ 
perature T(R) = 89 (R/10AU) -1 / 2 K, where R is the 
distance from the central star. This model is con¬ 
sistent with the distribution of the temperature and 
the surface density estimated from the observations of 
dust continuum (Beckwith et al. 1990). Hydrostatic 
equilibrium determines the density distribution per¬ 
pendicular to the midplane of the disk (e.g., Aikawa 
et al. 1996). 

We consider the gas-phase reactions, adsorption 
onto grains, and thermal desorption from grains. The 
reaction network we use is essentially the same as the 
one described in Aikawa et al. (1996). We adopt the 
UMIST94 data base (Millar et al. 1991; Farquhar 
& Millar 1993) for the gas-phase reactions, and also 
take into account some three-body reactions referring 
to Brasseur & Solomon (1986). Table 1 shows all the 
atoms and molecules included in our network (except 
H, H 2 , and He) and their adsorption energies. Be¬ 
cause the attenuation of cosmic rays is inefficient, we 
adopt the ionization rate in the dark cloud, ( « 10 -17 
s _1 . We neglect the ionization and dissociation by 
the interstellar and stellar ultraviolet radiation since 
it can be attenuated by grains in a thin surface layer 
of the disk (see § 4.3). We take the sticking proba¬ 
bility S = 0.3 for the collision of a gas particle with 
a grain (Williams 1993). For simplicity, we assume 
that all grains are spherical with radius a = 10~ 5 cm 
and that the dust-gas ratio is the same as that in the 
interstellar clouds. Although the dust-gas ratio near 
the midplane might be enhanced by sedimentation of 
dust, it does not affect much the results; we have ob¬ 
tained almost the same results with the dust-gas ratio 
higher by orders of magnitude. 
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Since our reaction network consists of one- and 
two-body reactions and contains only a small num¬ 
ber of three-body reactions, we consider the region 
of the density by number of hydrogen nuclei, «h <, 
10 12 cm -3 , which corresponds to i? ^ 10 AU for the 
minimum-mass solar nebula. With a typical rate coef¬ 
ficient k ~ l(U 30 cm 2 * * * 6 s _1 for the three-body reactions 
we have in this density range /c?z(H 2 )<^ 10 _ 18 cm 3 s _1 
even for the most abundant species H 2 as a third 
body. This is smaller than the characteristic rate 
coefficient ~ 10 “ 17 cm 3 s _1 for radiative association, 
which is the slowest among the two-body reactions. 
Thus the three-body reactions are not important in 
this density range. 

For simplicity and clarity, we assume that the den¬ 
sity and temperature do not change with time. Since 
we are mainly interested in the quiescent (WTTS) 
phase, we neglect the radial transport of matter and 
perform calculation up to ~ 10 7 yr. The initial molec¬ 
ular abundance is the same as in Aikawa et al. (1996), 
which assumed mostly the observed abundance in the 
dark cloud TMC-1. However, the molecular abun¬ 
dance may change during the active phase of the disk. 
To find out the influence of the initial condition, we 
also perform calculation with another set of the ini¬ 
tial abundance; all carbon is in the form of CO, the 
remaining oxygen is locked up in water ice mantle, 
and nitrogen is in N 2 . 

3. NUMERICAL RESULTS 

Since the result is found to be much more sensitive 
to the temperature than to the density, we show as 
the representative cases the results at T = 30 and 90K 
corresponding to the regions of R = 87 and 9.7AU 
with the density tih at the midplane 2.9 x 10 9 and 
1.2 x 10 12 cm -3 , respectively. 

Figures la and 16 show the results for the first 
set of the initial condition. The solid lines show the 
time variation of the abundance of carbon-bearing 
molecules relative to hydrogen. In the region of T = 
30K (Figure la) the abundance of CO remains nearly 
constant at its initial value n(CO)/nn ~ 7 x 10 -5 up 
to 10 6 yr. Thereafter CO decreases steeply, and CO 2 
ice and HCN ice become dominant. The main pro¬ 
cess of forming CO 2 and HCN is as follows. Through 
the reaction with H)j~, CO is transformed into HCO + , 
some of which become HCO mainly through the grain 
surface recombination. Reactions of HCO with C 


and N finally form CO 2 and HCN in the gas phascQ. 
Once CO 2 and HCN are formed, they are adsorbed 
by grains in a time as short as ~ 10 yr (Aikawa et al. 
1996). Since these molecules are hardly desorbed at 
T = 30K, they are locked up in the ice mantle, and 
as a result carbon is depleted from the gas phase. At 
somewhat higher temperatures (e.g., T = 90K as in 
Figure 16), the thermal desorption is efficient enough 
to compensate the adsorption; although CO 2 becomes 
abundant, CO is also abundant on a time scale of 10 7 
yr- 

The dotted lines in Figures la and 16 show the 
time variation of the abundance of oxygen-bearing 
molecules. Oxygen is initially assumed to be mostly 
in the form of O atom and CO molecule. The abun¬ 
dance of O 2 in the gas phase and H 2 0 ice in the 
grain mantle increases extensively with time. The 
H 2 O molecules form mainly in the following process. 
The reaction between O and Hg produces OH + , most 
of which are transformed into H 3 CU through H atom- 
abstractions from H 2 . Then HgO + recombines to 
produce H 2 0. In the later stages where O 2 becomes 
dominant, the main reaction of forming OH + is the 
H atom-abstraction by 0 + , which is produced via de¬ 
struction of 0 2 by He + . Since H 2 0 freezes onto grains 
at these temperatures while O 2 does not, most oxygen 
freezes out in the form of H 2 O ice and the abundance 
of O 2 decreases in a time scale of several x 10 6 yr. 

The dashed lines in Figure 1 show the time varia¬ 
tion of the abundance of nitrogen-bearing molecules. 
Initially most nitrogen is assumed to be in the form 
of N atoms. As time goes on, N2, HCN ice, and NH3 
ice become dominant. The main formation process 
of NH3 is as follows. The reaction of N2 with He + 
forms N + , which is transformed into NHj by repeat¬ 
ing H atom-abstraction. Finally NHj is turned into 
NH3 through dissociative recombination and grain- 
surface recombination. At temperatures below about 
70 K, NH3 freezes onto grains and most nitrogen is de¬ 
pleted from the gas phase. At higher temperatures as 
in Figure 16, NH3 cannot freeze and N2 is dominant 
even at later stages. 

4. DISCUSSION 


2 By recombination on grains HCO + is transformed into HCO 

radical or CO + H with unknown branching ratio. We as¬ 

sumed that the former branch is dominant. The time scale of 

the formation of CO 2 and HCN becomes longer if we assume 

otherwise. 
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4.1. Implications to Planetary Science 

We review the main points of our results and dis¬ 
cuss the implications to planetary science. First, we 
have found that CO and N 2 are transformed into C0 2 
and NH 3 . The main cause for this is in cosmic rays, 
which have not been taken into account in the pre¬ 
vious works. Cosmic rays produce chemically active 
ions such as Hg" and He + which can destroy CO and 
N 2 . This results in the depletion of CO in the gas 
phase and hence the faintness of its emission lines in 
the cold regions (T <; 70K) of the disks at the age ^ 
several xlO 6 yr. This also means that C0 2 ice and 
NH 3 ice will be detected in these regions. Secondly, 
our results show that oxidized (C0 2 ) ice and reduced 
(NH 3 and hydrocarbon) ice coexist in the protoplan¬ 
etary disks. This coexistence has also been suggested 
by the KI model (Prinn 1993). This result is to be 
noted because their coexistence is one of the most im¬ 
portant characteristics of comets (Yamamoto 1991). 
CO and N 2 would also coexist in ice mantle if they are 
physically trapped in water ice (Bar-Nun et al. 1985), 
although this effect is not included in our calculation. 
Finally, we have found that the molecular abundance 
in the gas and in ice mantle varies considerably with 
the distance from the central star. Comparison of 
this result with the molecular abundance of comets 
will enable us to deduce where comets formed, al¬ 
though quantitative discussion is out of the scope of 
this paper. 

4.2. Dependence on the Initial Condition 

The initial condition we have adopted so far is that 
considerable amount of carbon, oxygen, and nitro¬ 
gen are atomic (see § 2 in Aikawa et al. 1996). In 
order to check the dependence on the initial condi¬ 
tion, we have made calculation for another set of the 
initial abundance mentioned in § 2. Figure 2 shows 
the result for this case in the region of T = 30K. 
We can see that C0 2 ice becomes dominant in this 
case as well. Although there are no oxygen atoms 
initially, they are extracted from CO by the reaction 
with He + . The remaining carbon is transformed into 
hydrocarbon. At higher temperatures (T ;> 70K), 
C0 2 becomes abundant at later stages, although its 
abundance nco 2 /nH = 10 -5 — 10 -6 is lower than that 
for the first case on the initial condition. Similarly, 
N 2 is destroyed by He + and is transformed into HCN 
ice and NH 3 ice. 


4.3. Discussion on the Assumptions 

Although we have neglected the radial transport 
of matter, our results are qualitatively applicable to 
the disks with radial transport such as the accretion 
disk and weakly-turbulent disk, because the chemical 
paths described in § 3 are always active and domi¬ 
nant in the reaction network as long as the matter 
is in the region of low temperature and low density 
(R ;> 10AU in the minimum-mass solar nebula). The 
molecular evolution in the accretion disks will be de¬ 
scribed in the forthcoming paper in more detail. On 
the other hand, our results are not applicable to the 
disks with efficient mixing where the strong turbu¬ 
lence or convection mixes matter between the inner 
hot (T severalx 100K) region and the outer cold 
region (R ^ 10AU) in a time scale 10 6 yr, as inves¬ 
tigated by Prinn (1993). 

We have neglected the UV photolysis. Smaller 
grains sink more slowly toward the midplane of the 
disk (Nakagawa, Nakazawa, & Hayashi 1981). In 
addition, even a weak turbulence which cannot ef¬ 
ficiently transport matter radially would be sufficient 
to prevent the sedimentation of small grains. There¬ 
fore, it would be possible that small grains shield the 
main part of the disk from the UV radiation. Mid- 
to far-infrared excesses observed in the spectra of 
WTTS strongly suggest the existence of small grains 
at R ;> 10AU even in the relatively quiescent disks 
(Strom et al. 1989). 

If almost all grains sank to the midplane in less 
than about 10 6 yr despite the above discussion and 
the mid- to far-infrared excesses in WTTS, our re¬ 
sults were not applicable. The molecular abundance 
must then be affected by the photolysis because the 
line shielding is not efficient enough for some kinds of 
species. 
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like to acknowledge support by JSPS Research Fel¬ 
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Table 1 

THE ADSORPTION ENERGIES OF ATOMS AND MOLECULES IN KELVIN 


Species 

Adsorption 
Energy (K) 

Species 

Adsorption 
Energy (K) 

C 

800 

C 2 H 

1460 

N 

800 

c 3 

2010 

O 

800 

CCN 

2010 

Na 

11800 

co 2 

2690 b 

Mg 

5300 

NH 2 

860 

Si 

2700 

no 2 

2520 

S 

1100 

OCN 

2010 

Fe 

4200 

OCS 

3000 

HS 

1500 

o 2 h 

1510 

CH 

650 

o 3 

2520 

c 2 

1210 

so 2 

3460 b 

CN 

1510 

h 2 co 

1760 

CO 

960 b 

h 2 cs 

2250 

cs 

2000 

ch 3 

1160 

NH 

600 

c 2 h 2 

2400 c 

n 2 

710 c 

c 3 h 

2270 

NO 

1210 

c 3 n 

2720 

NS 

2000 

nh 3 

3080 b 

OH 

1260 

h 2 c 3 

2110 

o 2 

1210 

hc 3 n 

2970 

SiH 

2940 

ch 4 

o 

00 

o 

SiC 

3500 

ch 3 o 

1710 

SiO 

3500 

c 3 h 2 

2110 

SiS 

3800 

ch 3 cn 

2270 

so 

2000 

ch 3 oh 

4240 b 

h 2 o 

4820 b 

ch 3 o 2 

2370 

h 2 s 

1800 

c 3 h 3 

2220 

HCN 

4170 c 

h 3 c 3 n 

3270 

HCO 

1510 

c 3 h 4 

2470 

ch 2 

960 

CHsOOH 

2620 


a Allen & Robinson 1977 and Hasegawa & Herbst 
1993 except for 6 and c 

b Sandford & Allamandola 1993 
c Yamamoto, Nakagawa, & Fukui 1983 
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Fig. 1.— Evolution of the molecular abundance in 
two representative regions of the disk, (a) R = 87AU 
(n H = 2.9 x 10 9 cnr -3 , T = 30K) and (b) R = 9.7AU 
(riH = 1-2 x 10 12 cnr -3 , T = 90K). The initial molecu¬ 
lar abundance was determined referring to the abun¬ 
dance in dark clouds. The solid, dotted, and dashed 
lines represent the abundance of carbon-, oxygen-, 
and nitrogen-bearing molecules, respectively, relative 
to hydrogen. 

Fig. 2.— Evolution of the molecular abundance at 
R = 87AU (n H = 2.9 x 10 9 cnr" 3 , T = 30K) with the 
initial condition that carbon is in the form of CO, ni¬ 
trogen is in the form of N 2 , and the remaining oxygen 
is locked up in water ice mantle. The other details are 
the same as in Fig. 1 


REFERENCES 

Aikawa, Y, Miyama, S. M., Nakano, T., & Ume- 
bayashi, T. 1996, ApJ, 467, 684 

Allen, M., & Robinson, G. W. 1977, ApJ, 212, 396 

Bar-Nun, A., Hertman, G., Laufer, D., & Rappaport, 
M. L. 1985, ICARUS, 63, 317 

Beckwith, S. V. W., Sargent, A. I., Chini, R. S., & 
Glisten, R. 1990, AJ, 99, 924 

Brasseur, G & Solomon, S. 1986, in Aeronomy of the 
Middle Atmosphere, ed. Brasseur, G & Solomon, 
S. (D. Reidel Publishing Company) 

Farquhar, P. R, A., & Millar, T. J. 1993, CCP7 
Newsletter, 18, 6 

Hasegawa, T. I., & Herbst, E. 1993, MNRAS, 261, 83 

Hayashi, C. 1981, Prog. Theor. Phys. Suppl., 70, 35 

Millar, T. J., Rawlings, J. M. C., Bennett, A., Brown, 
P. D. & Charnley, S. B. 1991, A&ASS, 87, 585 

Nakagawa, Y., Nakazawa, K., & Hayashi, C. 1981, 
ICARUS, 45, 517 

O’dell, C.R., & Wen Zheng 1994, ApJ, 436, 194 

Prinn, R. G. 1993, in Protostars and Planets III, ed. 
Levy, E. H., & Lunine, J. I., (Tucson: Univ. of 
Arizona Press), 1005 

Sandford, S. A., & Allamandola, L. J. 1993, ApJ, 417, 
815 

Sargent, A. I., & Beckwith, S. V. W. 1994, Ap&SS, 
212, 181 

Strom, K. M., Strom, S. E., Edwards, S, Cabrit, S., 
and Skrutskie, M. F. 1989, AJ 97, 1451 

Umebayashi, T., & Nakano, T. 1981, PASJ, 33, 617 

Williams, D. A. 1993, in Dust and Chemistry in As¬ 
tronomy, ed. Millar, T. J., & Williams, D. A. (Lon¬ 
don: Institute of Physics Publishing), 143 

Yamamoto, T. 1991, in Comets in the Post-Halley 
Era, ed. Newburn, R. L. Jr. et al. (Kluwer Aca¬ 
demic Publishers), 361 

Yamamoto, T, Nakagawa, N., & Fukui, Y. 1983, 
A&A, 122, 171 

This 2-column preprint was prepared with the AAS iATgX 


6 



macros v4.0. 


7 












Cf 


>D!O z H 






